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<H conserved charges of the double copy model equations parametrising the correlators of twist 
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' ^ ' parametrising correlators of primary twist fields are compared to those of the literature and 
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1 Introduction 



Correlation functions of local fields of Quantum Field Theories contain all the physical in- 
formation of the model. These functions are often only accessible perturbatively and at 
large energies, but in two-dimensional, integrable quantum field theories with free fermion 
representations non-trivial correlation functions can be expressed as solutions to integrable 
differential equations. The first example of this phenomenon is the spin-spin correlation func- 
tion in the thermally perturbed Ising model [22] and other results, some concerning models 
with Dirac fermions, have followed through a variety of methods: holonomic quantum fields 
(Dirac model) [19], Fredholm determinants (Ising and Dirac models) [1, 11, 3], determinants 
of Dirac operators (Dirac model in fiat and curved space and with magnetic field) [15, 16, 12], 
and doubling of the model (Ising spin chain, and Ising QFT model at zero and non-zero tem- 
perature and in curved space) [17, 9, 6]. The existence of such differential equations allows 
the relevant correlation functions to be evaluated very accurately once initial conditions have 
been fixed via conformal perturbation theory and form factor analysis. 

Twist fields are local fields which have non-trivial monodromy, with respect to the fermion 
fields, associated to a symmetry of the model. From the holonomic quantum field descrip- 
tion we expect general correlation functions of twist fields to be parametrised by integrable 
differential equations. Twist fields were first introduced in [20], as the Z2 monodromy fields 
of the Majorana fermion corresponding to the spin field of the Ising model. Similarly, it 
is known that primary U{1) twist fields in the Dirac model reproduce correlation functions 
of exponential fields in the sine-Gordon model at a particular value of the coupling (the 
free-fermion point) [13]. It is correlation functions of such Z2 and U{1) primary twist fields 
that have been studied and that are known to lead to differential equations in QFT models. 

This paper is an extension of the work presented in [7] where the differential equations 
for the correlation function of twist fields with equal monodromy were found using a method 
similar to that of [9] . In [9] a double copy of Ising field theory is considered and its conserved 
charges used to derive Ward identities which lead to differential equations for the spin-spin 
two point function. The states considered in [7] were kept as general as possible and only 
translation and parity symmetry were required in order to reduce the Ward identities to the 
known form, [3]. In the present work we use the same methods to retrieve the differential 
equations parametrising the correlation functions of twist fields with differing monodromy. 
To obtain these equations it transpires that the states under consideration are also required 
to be rotation symmetric as without this property the Ward identities derived below do not 
provide sufficient information to obtain the result. 

The twist fields of the Dirac model, (t„ for a G M, are characterised by their U{1) mon- 
odromy, e^"^'" and dimension a^. The properties of these fields, along with their fermionic 
descendents, are discussed in more detail in [7]. The main aim of this work is to reproduce 
the differential equations parametrising the correlation function of twist fields. When written 
as 

(fT„(a;,y)a^(0,0)) = c„C;3m"'+^%^(-'^) (1) 

where E is a function of r = (1/4) (a;^ + y^) we find that it is the a solution of the following 
equations, fist found in [18]: 

^(l-cosh(2V^)) (2) 

^ sinh(2^) + tanh V(l - tanh^ V)- (3) 
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+ -dr ] V 
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In deriving these equations we find a discrepancy between (3) and the results of [3], 
namely a factor of 4 multiplying the (a — /3)^ is not present in our equations. As we are 
unable to find an error, or explain this missing factor, we attempt to verify which equation is 
correct. Using the form factor expansion of the correlator at large r we obtain an approximate 
expression for E and thus also for ip. We then examine how well this expression solves both 
(3) and the equation of [3]. Numerical and analytical methods are used to this end with 
each one indicating that the equations presented in this paper do indeed give the correct 
parametrisation. 

The method presented in this paper also gives differential equations for the correlation 
functions of descendent twist fields. 



(a„+i^„(a;,y)a^,^+i(0,0)) =ge'''(^-")c„c^m"'+"+'''+^+ie^'W (4) 

parametrised by the same auxiliary function ■0. Both correlators are shown to have the same 
r dependence and the constants p and q are discussed in more detail in section 4.1, but 
essentially depend on the positions of the fields through their braiding relations. It is shown 
the function E' satisfies 

sinh -0 cosh ■0 



r 



dt + -dr E' = ^ ~m'~ -\ ^ (5) 



which is a new result that appears naturally from the Ward identities. 

This paper is set out as follows: In section 2 notation is set out and the model under 
consideration is explained. We also present the fields and relevant form factors, however no 
derivation is presented and we refer the interested reader to [7] for more details. Section 
3 introduces the double model and the action of relevant conserved charges is described. 
Again more detail is presented in [7] . In section 4 the correlation functions of interest and 
their symmetries are discussed and the Ward identities derived from the conserved charges 
in the double model are written down. These Ward identities are then manipulated to give 
the differential equations parametrising the correlation functions of primary and descendent 
twist fields. As the equations for the primary twist fields differ from those of the literature, 
section 5 presents numeric and analytic evidence that the equations derived in this paper 
give the correct parametrisation. Finally section 6 gives our conclusions and outlook. 



2 The Model and its Fields 

In all that follows we will use the following conventions: a; S M is a space coordinate while 
y G M is the usual Euclidean time. These are combined into the complex coordinates z := 
— |(a; + iy), z := |(x — iy), with derivatives d = id^ + dy and d :— dg = —id^ + dy. 

When only one coordinate is given it is understood to be the space coordinate with time 
being set to 0. 

2.1 Dirac Fermions 

The standard, two dimensional, free Dirac Fermi field of mass m is an operator solution, 
y), ^'i(a;,y), to the equations of motion 

+ + 
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and the equal-time anti-commutation relations 

{■fR{xi),^''j^{x2)}^^7Td{xi-X2), {^'L(2:i),^'[(a;2)} = 47r(5(xi -X2) (7) 

with all other combinations anti-commuting. This is subject to the condition there exists a 
vacuum state, |vac), with the properties 

lim '^j^^L{x,y)\Ya.c) ~ lim ^I't, , (x, ?/)|vac) = 0. (8) 

The anti-commutation relations (7) have been chosen so as to give the fields the 'CFT' 
normalisation in terms of the z variable: as \zi — Z2I — ^ 0, 



(vac|r •^''j^{xi,yi)^B.{x2,y2) 
(vac|r [*[(a;i,yi)«'L(a;2,y2) 



I vac) ^ — — 



vac) 



Zl - Z2 

1 

Z1-Z2' 



Where T is time ordering, with operators at later Euclidean time being placed to the left of 
those at earlier times and picking up a sign, in general, if those operators are fermionic. 

These operators can be written in terms of creation and annihilation operators (Fourier 
modes): 

■^L{x,y) = J d9e-'^/^ {iDl{e)ey^'-'''P'' - D_{9)e-y'^''+'''P''^ (9) 

where 

Eg = m cosh 6*, pe = msinh^ (10) 
and their Hermitian conjugates are 



remembering that since y is Euclidean time it changes sign under Hermitian conjugation. 
The creation and annihilation operators satisfy the anti-commutation relations 

{Dl{ei),D+{e2)}^ 5(01-62), {Dl{9,),D^{e2)}^d{0i-02) (12) 

(with all other combinations anti-commuting), which are derived from (7), and the vacuum 
condition 

D±(6')|vac) =0 

from (8). The multi-particle states 

\0i---0„),,...,^:^ Dl{0i)---Dlj0r,)\Yac) for 0^ > . . . > e {+, (13) 

form a basis of the Hilbert space and have total energies Eg^ and total momenta pg. . 
Dual vectors are written as '^'^'"'^"{01 ■ ■ ■ 0n\ ■= \0i ■ ■ ■ fi'n)ei - £„ and clearly 

n 

<•< {0[--- 0'j0, ■ ■ ■ = n 4m - ^;)- (14) 



4 



It is useful to note that the identity operator can be written as 

^ 1 nOO poo 



2.2 Primary Twist Fields 

The internal C/(l) symmetry of the Dirac theory, "^r^l e^^^^^^fR^L, a G [0, 1), has associ- 
ated primary twist fields, aa{x,y). These are local, spin-less, U{1) neutral bosonic quantum 
fields with scaling dimension [19]. At this point we note that these twist fields can be iden- 
tified with exponentials of the sine-Gordon field at the free fermion point via ct^ =: e*"* :, 
where $ is the appropriately normalised sine-Gordon field. The relevant bosonisation rela- 
tions also give bilinears of fermions in terms of ^'J^^E'l =■ e** ■ and vE'i?^'{^ =: e~** :. 
Expressions for the fermions themselves require exponentials of terms which are non-local 
in terms of the sine-Gordon field and so the descendent fields discussed in section 2.3 are 
also given by non-local expressions, meaning that these fields are not so natural in the sine- 
Gordon theory. The implications of the identification with the sine-Gordon model are beyond 
the scope of this paper but discussed further in [3] . 

The twist property of CTq is manifested in the monodromy property of time ordered 
correlation functions 

Ciz) (vac|r[---*fl,,L(x,2/)a„(0)---]|vac). (16) 

When the complex variable z is continued along a loop 7, surrounding the origin once counter- 
clockwise, the function C is given by 

Cie^'^'z) = e-2'^*"C(z) (17) 

which holds for any loop which can be contracted to the origin without intersecting the 
position of any other twist fields present in the correlation function. When "^r^l is replaced 
with ^'^ ^ a similar property holds but with e~^'^*" replaced by e^'^'". This twist property 
can also be expressed via the exchange relations 

f fT„(0)*fl.L(x) (x<0) 

^rM^MO) - { 2.,^ _ (18a) 



e 



^a(0)*fl,L(2:) {X > 0) 



r cr a{0)'Sl Ax) (x<0) 
¥^r{x)a^iO) = { ''''' ' (18b) 

I e-'--a^ml,{.) [x > 0). 



The Hermitian conjugates of the primary twist fields is given by 



(19) 



From (18) the form factors of the primary twist fields can be calculated [10, 14, 8]. The 
fields are U{\) neutral and so only have non- vanishing form factors with U{1) neutral states: 

(vac|a„(O)|0i02...M+...+_..._ =c„m"^(-l)"("-i)/2 ('^^)" 

X {i{{u.)'^^-{u.,.)^-A "'<^^;^r"^]"T"'<^^7""^^ (20) 
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where there are n — 's and n +'s in the state, and Ui := exp(0i). In particular, the two-particle 
form factor is 

(vac|a„(O)|0i02)+- = c„m" ^ . ^ ^ ■ (21) 

Sttz cosh ^ 

These form factors may differ from those written in other publications in the choice of sign 
of a and the two particle form factor. The constant Cq, does not influence the calculations of 
this paper and so we direct the interested reader to [2], [13], [4] and [7] for further discussion 
and evaluation of this constant. 

Other matrix elements can be obtained by analytic continuation in the rapidities, as is 
standard in the context of l+l-dimensional QFT [21]. Note that the form factors (20) and 
the Hermiticity relation (19) are in agreement with the analytic-continuation formula 

(vac|a„(O)|(0i +nT)--- (02„ + i7r))+...+_..._ = +-+-- {e2n ■ ■ ■ ei|a„(0)|vac), (22) 

where the analytic continuation is simultaneous on all rapidities. 



2.3 Descendent Twist Fields 



As well as the primary twist field described in the previous subsection our model also contains 
descendent twist fields. The two families of descendent twist field of interest here occur in 
the operator product expansions of primary twist field with the Dirac fields 5'^^ and 4'!^: 



r 



(-iz)-"a„_i,„(0). 



(23) 



These OPEs are valid for all a E M\Z* (where Z* := Z\{0}) and the same families of 
descendent twist fields occur in the OPEs of ^I^^ and 'I'^ with primary twist fields: 



T[^l{x,y)aa{0) 



-(iz)-V„,„_i(0) 

- {iz)°' a a, a+l{0)- 



(24) 



The descendent twist fields cra±i^a have dimensions a'^ ±a+l/2, spins ±a-|-l/2 and charges 
=F1, as discussed in [7]. The one particle form factors of these fields can be evaluated from 
the OPEs: 



(vac|cra+i,a(0)|6') + 
(vac|f7Q_i,a(0)|6l)_ 



?7ra/2 

'"r(l + a) 



^Q2+a+l/2g(Q+l/2)S 



-ICo 



-a+l/2g(-Q-|-l/2)e^ 



r(i-a) 

All higher particle form factors can be calculated via Wick's theorem. 



(25) 
(26) 



3 The Double Model 

As was the case in [7] inspiration is taken from [9] where a model containing two non- 
interacting copies of Ising field theory is studied. The conserved charges in this double 
model allows differential equations for the Ising spin-spin correlation function to be written 
down. 
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Here we consider a double model consisting of two non-interacting copies of the Dirac 
model discussed in section 2 with and $ denoting the fermion fields and and e'^^ 
the respective creation and annihilation operators of the two copies. As these two copies do 
not interact their fermionic fields anti-commute. The copy which other fields belong to will 
be identified by a superscript; so twist fields from the two copies are denoted ct* and u*. 



3.1 Conserved charges 

Energy and momentum are conserved charges of the single-copy model. They are associated 
to the dynamical invariance under space and time translation: the equations of motion (6) 
possesses this invariance. In the double-copy model, the energy and momentum conserved 
charges are the sums of the corresponding charges of each copy. However, there are other 
conserved charges in the double-copy model that can be constructed out of these. As the two 
copies are non-interacting, the energy-momentum operators for each are still independently 
conserved quantities, and can be combined differently to give new conserved quantities. The 
two specific conserved charges that are of interest to us are the differences of those of the 
single copies. We define P and P via the following action (these conserved charges are chosen 
to be anti-Hermitian) : 

[p, 0*0*] =iao*o* - io*ao* 

[P,0'^0'^]=-idO'*0'^ + iO'^dO'^. (27) 

This holds for any local fields O* and O* interacting non-trivially with fields in copy ^E* and 
$ respectively. In particular, the actions on the creation and annihilation operators are 

[P,D±{0)] = -ime^D±{0), [P, E±{e)] = ime'^E±{e), 
[P,D±{e)] = ime-^D±{e), [P,E±{e)] = -ime-^E±{e) 

from which it is simple to derive an explicit expression for P and P through bilinears in the 
creation and annihilation operators. 

In the double-copy model, there is another conserved charge Z, related to the 0(2) 
rotation symmetry amongst the copies. Written in terms of the Fermi fields ^' and $, the 
charge Z is (again, chosen to be anti-Hermitian) 

^ = ^ / M'i'R'^R + + *^^$J^ + *i*L). (28) 

The action of this charge on the creation and annihilation operators is 

[Z, Di{6)] = -Ei{0) [Z, Ei{6)]=Di{e) (29) 

and similarly for D±{6) and E±{9). 

Higher order conserved charges can of course be obtained by commutations of P, P and Z 
but as was the case in [7] the Ward identities coming from these three charges prove sufficient 
to derive the equations of interest. Below, we will derive the integrable differential equations 
from the Ward identities associated to Z, and use the simpler space-time dependence of 
correlation functions coming from translation, parity and rotation symmetries to achieve our 
result. 

As was observed in [5] (in the generalised situation of a theory on the Poincare disk), it 
is the relation 

[P, [P, Z]] = [P, [P, Z]] = Am^Z (30) 
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which provides the dependence on the mass in the differential equations. Relation (30) 
corresponds essentially to the equations of motion of the theory, as expressed using the 
charge Z. It can be derived from (27) and (28) using the equations of motion (6). Relation 
(30) is the only one where the massive theory is used: besides it, we only need the action of 
Z on products of twist fields and their first derivatives in order to find differential equations 
parametrising the correlation function (1). 



3.2 Equations for the action of Z 

The action of the conserved charge Z on products of local fields 0*0* gives linear com- 
binations of similar products of local fields (with unchanged locality index) whenever the 
combined locality index of the product 0*0* is one. Otherwise, the result of the action is 
non-local, and the resulting Ward identities are not useful. The equations that are needed 
here are derived in [7] so we simply report the necessary results. 

In the rest of this section we neglect to write the location of the fields, to make the 
equations easier to read, it being understood that all fields are evaluated at the same point, 
(0, 0) for instance. The equations representing the action of the charge Z on twist fields 
involving no derivatives are: 





= 






■ i ^ 






■ I ^ $ 


_* ^ 


^* 1 


./ \]> $ \]/ 




derivative 


are as follows: 






«((5cr*„+i)cr*+i_„ 


- 0'a + l,a(5c^*a+l)) 


1^* 1 


*('^a,a+l(^cr*+l,Q) 


- (^cr*+l,Q)c^a,a-Hl) 






~o-*_i^„(5o-*„_i)) 




i(o-*a-l(^crLl,a) 


- (^CT*_i,„)cr*„_i) 



(<9cra+l,a)'^*a+l] = ^^^l^^^t+l) ^ {d<^l+\)ot) 

(^cr*+i,„)cr*„+l] = i((9cr*)f7*+i -CT*+i(9f7*)). 

The only relations involving second derivatives of twist fields that we need can be obtained 
by using those above along with the equation of motion of the charge Z (30). Indeed, we 
only need action of Z on double-derivative fields of the form [P, [P,0*0*]], which we can 
evaluate using 

[Z, [F,[P,0*0*]]] 

= [[P, [P, Z]] , o*o*] + [P, [Z, [P, 0*0*]]] + [P, [Z, [P, 0*0*]]] - [P, [P, [Z, o*o*]]] 
= 4to2[z^O*0*] + [p^ [P,0*0*]]] + [P, [Z, [P,0*0*]]] - [P, [P, [Z,0*0*]]]. 

This is how the mass dependence will appear in the equations for correlation functions. 
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4 Correlation Functions 

The aim of this section is to generahse the results of [7] and obtain a pair of partial differential 
equations parametrising the correlation function 

(a„(x,y)a^(0,0)) (31) 

using the Ward identities derived from the action of Z. 

Before describing the properties of the correlation functions we first introduce the notation 

K'fi^x,V) = {ao.,p{x,y)a^,,{^M (32) 

where it is understood that when the two indices are equal the twist field is a primary twist 
field, i.e. 



4.1 Symmetries of Correlation Functions 

Before writing down our Ward identities we will first discuss some symmetries of the cor- 
relation functions involved. Firstly since our system has translation symmetry it is always 
possible to place one of the fields at the origin, as in (32). 
Next we note that 

FtX\tX\{x.v)^Fti{x.v) (33) 

which follows directly from the definition of the primary twist fields and can be seen explicitly 
by inserting the identity (15) between the fields on each side and using the form factors (20). 
We also note that as the fields involved in F^'^ are spin-less we do not expect to this function 
to depend on the angle 9, where z — re** and z — re~**, so we may also write 

Fi:i{x,v)^Fii[r). (34) 

Along with translation symmetry we also assume parity symmetry, discussed in more 
detail in [7]. Combining this with the fact that the primary twist fields have bosonic statistics, 
we see that F^'^ (a;, y) — F'^''^ (a;, y). In order to achieve a similar relation for the descendent 
twist fields we need to return to the braiding relations (18) and recall that by definition 

o-a+i,a(fi') = lim (-i(z - w)y^^\{z)aa,ot(w) 
o-Q,Q+i(w) = \m\{i{z-w)Y°'^i^{z)ao,,a[w) (35) 

where the factors (^—i{z — i(j))~" and {i(z — iD))^" are taken on the principal branch and so 
are continuous exactly where ^'^(z)(Tc[,q(u>) and Liz)^ a,akw) are continuous. Using these 
definitions and the braiding relations we obtain 

(cr/3,/3+l(w2)cra+l,a(wi)) = <^ {wi)(Jp^p+i{w2)). (36) 

We can now write 

{aa+iAwi)crf3,fi+i{w2)) {i{wi - W2))'^~"/f+iti(|wi - W2\) (37) 
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which reproduces the correct rotation properties. The choice of the factor {i{wi — i«2))^~" 
is in some sense arbitrary, but defined on the principal branch, it is convenient as we can 
also write 

{ap^p+i{w,)a^+,,Uw,)) - -e2-C(j(«;i _ W2)f~v!^:,\+l{\w^ - W2\) (38) 

where the same function — W2\) is involved and C = a if x{wi) < x{w2) and C = P 

if x{wi) > x{w2)- This choice will not be a problem as we will always consider products of 
the form F^tt+i^a+i'a ^^'^ so the same constant will occur in every term. 
With all this in mind we can now write 

F^^f^}{x^v) = qe^'^f'-'^)h{T) (39) 

where p and q are phases, as in (38), coming from (18), and the functions /i and /2 will be 
discussed in the next subsection. 

4.2 Ward Identities 

We are now in a position to write down the Ward identities associated with the actions of 
P, P and Z and simplify the resulting equations. 
The Ward identities of interest here are: 

([Z,f7*Jz)a*+i,„+i(z)a*+i,^(0)a*^+i(0)]) = (40a) 

([Z, [P,<jz)a*+i.„+i(z)]a*+i,^(0)a*^+i(0)]) = (40b) 

([Z, Jz)a*+i,„+i(z)]a*+i,^(0)a|,^+i(0)]) = (40c) 

([Z, [F, [P,<jz)a*+i,„+i(z)]]a|+i,^(0)a|,^+i(0)]) = (40d) 

([Z, [P,<Jz)a*+i,„+i(z)][P,a|+i,^(0)a|,^+i(0)]]) = 0. (40e) 

Maintaining the full (z, z) dependence for the moment these identities become, respectively, 

F^XYMF'J.^li-) - F^^X^'^^\z)F^J,^^^M + F^i{z)Ff:XliXli-) = « (41a) 

dF'^%Y,{z)F^Jtli^) - F^^%Y.i^)dF^JtJz)+ 

dF^X'''^\z)F^J,^^^,{z) - Pfi^^^+^(-)5Ff;l.+i(.) = (41b) 

5<^„^?(.)Ff;^^(.) - F^^%Y.i^)-dF^Jtli^) 

- ~dF^X'^^+\z)F^J,^^^,{z) + F^^X'''^\z)-dF^J,^^_,,{z) = (41c) 

^^F^.lY^i■AF^t:'.{z) - <t^?(.)apf;^,^^(z) - BF^^lY.{^)dFt^.:i{z) 

+ F'^%YMdBF^Jtl{z) - ^m-F^^%YMF^ttl(z) + ddF^t^^^^\z)F^:,^^^^,iz) 
- dF^^X^-f>^\z)-dF^J,^^^,{z) - dF^X'''^Hz)dF^:l^,^^^,iz) + <r^^+^(^)9aFffi,„+i(z) 

- ddF^:^iz)F^tlitli-) + dF^£i^)dF^tlitli-) + ~dFt,^{z)dF^ltliXl{z) 

- ;f {z)ddF^„Xlill (z) = (41d) 
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+ F^%YMddF^J,^^l{z) + ddF^;^iz)F^t'dtl(-) ~ dFiiiz)dF^tlitli-) 
- dF!;^iz)dF^:iitli^) + F^:^{z)ddF^tlitli-) + ddF^%'^^^\z)F^S^^^M 
- dFP%^'^+\z)-dF^S^^^,{z) - dF^^X'^P+\z)dF^S^^^,{z) 

+ F^X'^^+\z)ddF^J,^^^,{z) = 0. (41e) 

Now switching to polar coordinates and inserting the expressions (39) into (41b) and (41c) 
we find that /i = /2 up to an irrelevant constant. Thus we let 

/(r):=/i(r)-/2(r). (42) 

We now depart from the method of [7] ; in the present case we have no information about 
the relation between F^^'^'^^ and F^]^^ ^j^-^ so we will use (41a, 41b, 41c) to eliminate these 
functions from (41d) and (41c). In doing this we will also simplify our notation, using / as 
in (42) and (39), and letting 

Ft,^{r)^F{r). (43) 

It is also simpler to switch to r and 9 coordinates, as the 9 dependence will drop out, and 
introduce the operator 

V^dl + -dr. (44) 
r 

Putting all this together we see that (41d) becomes 
- pqfVf + pg ("^/)' /2 + pq^2j2 ^ ^q^pY 

- p2 Ipqp {F'idrFf - 2pqFfdrFdrf + {pqff{drff - (pqf (a - Pfflr^) = 0. 

(45a) 

and similarly (41e) becomes 



-r pq 

1 



- pqfVf + pq'—f^f + FVF 

{F\drFf - 2pqFfdrFdrf + {pqf f (Or ff - {pqf {a - /Sff^) = 0- 



F^-pqP 



(45b) 



Now as in the case for a = (3 we can make progress by setting 

F + \fpq! — cosh((^) F — ^/pqj ~ sinh(cp) (46) 

which reduces (45a) and (45b), respectively, to: 

{Vx - X'<y9)(sinh(2(p) - cosh(2</?)) 

+ O.V^)^(2sinh(2^)^ (^"-^^;^ + ^^f 

cosh (p smh ip 

+ m^(cosh(2(/3) - smh{2ip)) 
+ ^^"p (cosh(2v?) - sinh(2(^))(l + ^(coth(p - 1)(1 - tanh^)) = (47a) 
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sinh(2(^)2?X + cosh{2(p)'Dip 

+ {dr(pfismh{2tp) - cosh(2(p) coth(2(^)) 



Eliminating x from these equations leaves 



(cosh(2(^) - sinh(2<y9))(coth(2(p) + 1) = 0. (47b) 



Vip - idrip)^{l + cotli(2v3)) - sinh(2(p)(sinli(2(p) - cosli(2<^)) 

4^2 



(l + coth(2(p)) =0 (48) 



which we attempt to transform into (3) meaning we must eliminate the dr^p term. This is 
accomplished by noting that 

h" 

— = -1 - coth(2v?) (49) 

is solved by 

2^ = M^)-ln(^5_^j (50) 
and it is then a straight forward exercise to show that 

™2 / _ o\2 

Vi,^— sinh(2V') + ^ tanh(V') (l - t&nh^{il))) . (51) 

To obtain (2) we note that by setting 

= 2F(r) = e^+'^ (52) 

we find 

r»S = r»x + '^'P 

= — (1 - cosh(2V')) . (53) 

As noted in the introduction (51) disagrees with the equations in [3] by a factor of 4 in 
the (a — /3)^ term, which is not just a choice of coordinates since m scales like 1/r. Evidence 
that (51) gives the correct parametrisation is presented in the next section. 

It is also worth noting that by setting 

e^' = 2^/(r) = e^-^ (54) 
we find that / is also parametrised by the same function, -0, via the equation 

V^'^^^-m^- (55) 
sinh ^ cosh 

which, as stated in the introduction, is a new result that appears naturally when applying 
this method. 
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5 Analysis of Differential Equations 

In order to show that the results of the previous section give the correct parametrisation of 
the correlation functions we examine the form factor expansion of the correlator 

(a„(x,y)a^(0,0)) (56) 

at large distances and look at how well this solves the differential equations (51). It turns 
out to be sufficient to only consider terms up to two particles so we use the expression for 
the two particle form factor given in (21). 

To begin we insert the resolution of the identity (15) in between the fields of (56) and 
then expand the two particle term: 

{a a{x,y)af3 (0,0)) = (vac|crQ,(a;, y)|vac) (vac|cr^(0, 0)|vac) 

d(?id02 (vac|a„ (x, y)\ei, 62) (^2, Oi\c7p (0, 0) | vac) + • • • 

(vac|CTQ, {x, y) I vac) (vacjcr/j (0) | vac) 

X ( (vac|cra(0,0)|6»i,6l2)+_(vac|a/3(0, 0)16*2 + i7r,di + iTr)+- ) + ■■■ 
CaCfjm" '^^ sin(7ra) sin(7r/3) 



47r2 



g(ei-82)(a-/3) 



X / d6'id6i2e™^(^''i+^''2)-^^(p«i+p«2) 2-e73e^ + ' ' ' (^7) 

\ .7 cosh ( 2 ) / 

Since the first term is constant it is the second term which will eventually provide the leading 
large distance behaviour of ip. First, it is useful to use boost invariance to replace the {x,y) 
dependence in this term with r dependence, so that the second term in (57) becomes 



c„c^™"+^ sin(7ra)sin(7r/3) /■ e(«-«^)(«-« 



4^2 J -'-'^ cosh2(^) 



(58) 



This expression can be further simplified using a change of variables inside the integral. 
Using the variables 



1 "^2 , , (71 - C2 , , 

1 = — ^ — and 02 = — - — (59) 



2 " 2 

we see that the integral in (58) becomes 

1 r p202(a-^) 

r\rh p-rm{co8h{cj>i+<l>2)+co8h{4>i-cl>2)^ 

L<-i(f'2t; 2 / 

cosh 02 

1 r p202(a-^) 



, rlr/)nP~2rm cosh((/)i) cosh(02) 

^ J cosh 



2 

g202(a-/3) 

2 2 i^o(2?'mcosh02) (60) 

cosh 02 



where Kn{x) is a modified Bessel function. 
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Armed with this expression we wish to find out what S looks hke at large values of r. 
From the definition (1) we can write 

, , sin(7ra) sin(7r;3) /• e^^C^-^ ^ ^^^^ 

E ^ In 1 ^ ' ^ ' / AO ^iCo 2mrcosh0 + • • • 

47r^ j cosh 9 

= HI - g{r)) + ■ ■ ■ (61) 
where g{r) is defined implicitly. For ease of notation we also introduce the constant 

sin(7rQ:) sin(7r/3) 



^ 47r2 
To approximate we need to calculate 



(62) 



dl +^dr]^ (63) 



r 



which quickly becomes very cumbersome when using (61). However, at large distances g{r) 
becomes small so we may expect the Taylor series expansion of the log to still provide a good 
approximation. Thus we define our approximate solution for S as 



-g{r) (64) 



and thus 



-dAt = -{g"{r) + -g'{T)) 
r J r 



4m2^ / de e^'^'^°'-'^^Ka{2mrcos\ie). (65) 



We now have a useful approximation to the left hand side of (53) but to find a usable 
expression for ^ we must also examine the right hand side. As r — )■ cxi (56) is expected to 
tend to a constant, namely CaCpm"' '^^ , and so in this limit it must be that V' ~^ and we 
may Taylor expand the right hand side of (53) and only keep the first term: 

2 

^ (1 - cosh(2V')) = -mV^ H (66) 

Letting our approximate solution be ?/; we see that 

^2^4^ I de e^''^°'-^^Ko{2mrcoshe) = J{r) (67) 



where the function J{r) is defined for later convenience. Our goal now is to see how well this 
approximate solution solves the equation (51). For comparison we will also examine how 
well the corresponding equation from [3] , 

= ^ sinh(2V') + "^^""^^^ tanh(?A) (l - ia.nh^{iP)) (68) 

is satisfied V'- 
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5.1 Numerical Approximation 



To give a first indication as to which solution is correct we will examine how well the ap- 
proximate solution (67) satisfies the equations (51) and (68) numerically. As the left hand 
side of both equations is the same it will be sufficient to examine the functions 

[dl + -dr)i) (69) 
r 

— sinh(2V^) + — ^tanhV'(l-tanh2 7/;) (70) 

— sinh(2?/;) + tanh?/i(l - tanh^ v3). (71) 
2 

These functions can be calculated numerically for a fixed value of {a — f3) and their curves 
plotted. The value of (a — /?) must be non-zero as the two equations are equal for a = l3. 
In figure 1 these curves are plotted with (a — /?) = 0.3 and r G (7, 7.25). From this plot it 
is clear that for modest values of r the expression presented in section 4.2, (71), is a better 
approximation of (69) than (70). The same behaviour is observed for different values of 
(a-/3). 



0.0001 05 r 




Figure 1: Plot of (69) (blue), (70) (yellow) and (71) (red). 

From plots such as figure 1 we see that (71) provides a better approximate solution for 
all values of r. At r = 10, for example (69) is 4.36163 x 10"^, (71) is 4.34464 x IQ-^ while 
(70) is 4.30542 x 10^^. So in this case and using the equation presented in this paper the 
approximate solution ip is out by 0.39%. When the equation (70) is used this difference 
jumps to 1.39%. Again the same pattern appears when different values of (a — /?) and r are 
used. 

The numerical evidence presented in this subsection strongly suggests that the equation 
presented in this paper is indeed the correct parametrisation for the correlation functions 
twist fields. 
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5.2 Analytic Approximation 

With the strong numerical evidence pointing to (51) being the correct expression for the 
function parametrising the correlation functions of twist fields we will now look for some 
analytic evidence to support this. The objective is to expand both sides of (51) and test 
whether or not the leading r terms cancel when the approximation ip = ip is applied. 

To begin with, as we know that "0 is small at large values of r, we may expand the right 
hand side of (51) and only keep the linear term: 



„ sinh(2?A) + ^ tanhi/'(l - tanh-" ?A) = + ^ V H (72) 

2 V / 

We now want to insert our approximation of ip and examine the equation 

{dr + -dr)^ - + + • • • (73) 



It is convenient at this point to use the function J{r) defined in (67) so that this equation 
becomes, neglecting higher order terms, 

J"J-\{J'f + lJ'J ^ ^ {a^(3f 



2J3/2 + (74) 

and to further ease our calculations this expression can be rearranged to give 

J" J - -( J')" + -J' J - 2 fm^ + 'f^' ) J\ (75) 
2 r \ T'^ ) 

In order to progress further we need: 

J(r) = 4^ J de e^'^^°'-'^^Ko{2mrcoshe) (76) 

J'(r) = -8m^ J de e^^"-^^^ cosh9Ki{2mrcoshe) (77) 

J"{r) = l&rn^i j d6 eA"''^^'^ cosh^ eKo{2mrcoshe) 

+8me j de e2("-'5)«^^Xi(2mrcosh0). (78) 

From this we see that the second term in the J" J term of (75) cancels with the J' J /r term 
so there are only two terms left to deal with. 

Next we note that there are no {a — P)"^ terms in these derivatives so we need to eliminate 
these from (75). This is done by integrating J{r) by parts, integrating the exponential term 
to bring down l/(a — /3) terms: 

J de e^^°''^^^Ko{2mrcosh9) = --^^^ J dO e^'-'^^'^'^^Ki{2mrcosh9)smh9 

d9 e^^""-^)" {2mrKoi2mrcoshe) sinh^ 9 



mr 



2{a- PY 
) 

- Ki {2mr cosh 9) — Ki {2mr cosh 9) cosh 9 



sinh^ 9 
cosh^ 



(«-/3)- 



j d9 e^^"-^^'^ Ko{2mr cos\\9){cosh^ 9 - I) 



^d6' e2("-'3)«^'— i^i(2mrcosh6'). (79) 



2(a-/3)2 J cosh 6* 
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Now writing the right hand side of (75) as 



^2j2 ^ J2 (80) 

and using the above integration by parts to rewrite one of the J's multiplying the (a — /3)^ 
term we see that the term is cancelled by part of the first integral in (79) and the other 
part of this integral cancels some of the integral from the J" J term on the left hand side. 
So bringing all the remaining terms of (75) together we have 



X ( de e2(""^)''( 5— — ifi(2TOrcosh6l) +cosh^6'is:o(2TOrcosh6l)) 

\ J 2mr cosh & 

-32to2^2 (^J de e^^""'^)^ cosh 6lifi(2mr cosh ( 



e)j c^O (81) 

from which we want to extract the leading large r behaviour. Expanding the Bessel functions 
for large r gives leading terms proportional to l/\/r: 



^ — 2mr cosh 6 



Amr cosh 9 



X ( / d0e2("-««cosh2 0j e-^"'-™^^^ 

4mr cosh 9 



-Z2m^e{ /dge"("-^)^coshg./ ^ -2m.cosh9\ (82) 
\J V 4TOr cosh6' / 

While these terms do not cancel directly we observe that the main contribution from each 
integrand is from the region around = and so the first order saddle point approximation 
of these integrals does indeed vanish, as we had hoped. 

6 Summary 

As stated in the introduction we have derived the equations parametrising the correlation 
function 

(a„(x,y)a^(0,0)) (83) 

following the method of [7]. We have also presented strong evidence that the equations 
derived in this paper provide the correct parametrisation of the correlation function. 

With this derivation we have retrieved the known results regarding twist fields in the Dirac 
model using a direct and straight forward method. Unlike the equations of [7] the equations 
found here require rotation symmetry and so the result is not as strong. One extension of 
this work would be to determine if it would be possible to derive similar differential equations 
where the states considered have more general symmetries. 

The method presented has already been applied to Ising field theory [9] but its success 
in this more general setting suggests that it may have wider applications. It would be 
interesting to discover if this method could be successfully applied to other massive field 
theories to produce new results. 



17 



Acknowledgements The author wishes to thank Benjamin Doyon for his continued 

support and guidance and Ed Corrigan and Peter Bowcock for their help in preparing this 
manuscript. This work was conducted under an EPSRC studentship. 

References 

[1] Babelon, O. and Bernard, D.: From form-factors to correlation functions: The Ising 
model, Phys. Lett. B 288 (1992) 113 

[2] Basor, E. L., and Tracy, C. A.: Asymptotics of a tau function and Teplitz determinants 
with singular generating functions. Int. J. Mod. Phys. A 7S1A (1992) 83. 

[3] Bernard, D. and LeClair, A.: Differential Equations for Sine-Gordon Correlation Func- 
tions at the Free Fermion Point, Nucl. Phys. B 426 (1994) 534; Erratum ibid. B 498 
(1997) 619 

[4] Doyon, B.: Two-point functions of scaling fields in the Dirac theory on the Poincare 
disk, Nud. Phys. B 675 (2003) 607 

[5] Doyon, B.: Correlation Functions in Integrable Quantum Field Theory, Ph.D. thesis, 
Rutgers University (2004). 

[6] Doyon, B. and Fonseca, P.: Ising field theory on a pseudosphere, J. Stat. Mech. (2004) 
P07002 

[7] Doyon, B. and Silk, J.: Correlation functions of twist fields from Ward identities in the 
massive Dirac theory, J. Phys. A 44 (2011) 295402 

[8] Efthimiou, C. and LeClair, A.: Particle - Field Duality And Form-Factors From Vertex 

Operators, Commun. Math. Phys. 171 (1995) 531 

[9] Fonseca, P. and Zamolodchikov, A.: Ward identities and integrable differential equations 
in the Ising field theory, arXiv:0309228 (2003) 

[10] Karowski, M. and Weisz, P.: Exact Form-Factors In (l+l)-Dimensional Field Theoretic 
Models With Soliton Behavior, Nud. Phys. B 139 (1978) 455 

[11] Korepin, V.E., Bogoliubov, N.M. and Izergin, A.G.: Quantum Inverse Scattering 
Method and Correlation Functions, Cambridge Univ. Press (1993) 

[12] Lisovyy, O.: On Painleve VI transcendents related to the Dirac operator on the hyper- 
boHc disk, J. Math. Phys. 49 (2008) 093507 

[13] Lukyanov, S. and Zamolodchikov, A.: Exact expectation values of local fields in quan- 
tum sine-Gordon model, Nucl. Phys. B 493 (1997) 571 

[14] Marino, E. C, Schroer, B. and Swieca, J. A.: Euclidean Functional Integral Approach 

For Disorder Variables And Kinks, Nud. Phys. B 200 (1982) 473 

[15] Palmer, J.: Determinants of Cauchy-Riemann operators as r-functions. Acta Applican- 
dae Mathematicae 18 (1990) 199 

[16] Palmer, J., Beatty, M. and Tracy, C. A.: Tau functions for the Dirac operator on the 
poincare disk, Commun. Math. Phys. 165 (1994) 97 

[17] Perk, J. H. H.: Equations of motion for the transverse correlations of the one- 
dimensional XY-model at finite temperature, Phys. Rev. A 79 (1980) 1 

[18] Sato, M., Miwa, T. and Jimbo, M.: HOLONOMIC QUANTUM FIELDS. 3, (1978) 

[19] Sato, M., Miwa, T. and Jimbo, M.: Holonomic Quantum Fields. 4, Publ. Res. Inst. 
Math. Sci. Kyoto 15 (1979) 871 



18 



[20] Schroer, B. and Truong, T. T.: The order/disorder quantum field operators associated 
with the two- dimensional Ising model in the continuum limit, Nucl. Phys. B 144 (1978) 
80 

[21] Smirnov, F. A.: Form-factors in completely integrable models of quantum field theory, 
Adv. Ser. Math. Phys. 14 (1992) 

[22] Wu, T. T., McCoy, B. M., Tracy, C. A., and Barouch, E.: Spin spin correlation functions 
for the two-dimensional Ising model: Exact theory in the scaling region, Phys. Rev. B 
13 (1976) 316 



19 



